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Abstract 

o . . 

In inflation cosmologies, cosmic structure develops through the gravitational instabil- 
^Dity of the inevitable quantum noise in primordial scalar flelds. I show how the acceleration of 
the universe deflnes the shape of the primordial spectrum of gravitational metric and scalar 
i^fleld fluctuations. I assess how we can determine the shape and overall amplitude over the 
i^flve decades or so of spatial wavelengths we can probe, and use current data from cosmic 
ground radiation (CMB) anisotropies, large scale clustering and streaming observations, 
^distribution functions of cosmic objects, to show how far we are in this program. Broad-band 
^ power amplitudes are given for CMB anisotropy detections up to spring 1994, covering angular 
C/^ scales from all-sky down to arcminutes: DMR, FIRS, Tenerife, SP91, BigPlate, Python, SP89, 
■^MAX, MSAM, ARGO, White Dish, OVRO. It may be only a little premature to say that a 
spectral shape is emerging which is near that of preferred inflation models. The cluster-scale 
density fluctuation power derived from the broad-band power for COBE must fall within a 
narrow range to get the abundance of clusters right. This rules out many structure formation 
models, in particular restricting the primordial spectral index to be close to the scale invariant 
one. I show that COBE band-powers found with full Bayesian analysis of the 53, 90, 31 a+h 
GHz flrst year DMR (and FIRS) maps are in good agreement, and are essentially independent 
of spectral slope and degree of (sharp) signal-to-noise flltering. Further, after (smooth) opti- 
mal signal-to-noise flltering («, e,, Weiner-flltering), the different DMR maps reveal the same 
large scale features and correlation functions with little dependence upon slope. However, the 
most probable slope depends upon how the maps are flltered: with no flltering whatsoever, 
the slope is high, but the power is not described by a single-slope law; as flltering is increased, 
the index moves nearer to inflation predictions.^ 



1 in Relativistic Cosmology, ed. M. Sasaki, Academic Press (1994); Proc. 8th Nishinomiya- Yukawa Memorial 
Symposium, Japan, October 1993 
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1. Introduction 



The inflation paradigm — that the region of the Universe in which we reside was once 
in a state of accelerated expansion — remains the best way to account for local homogene- 
ity and isotropy. And, of course, quantum noise generated during acceleration is a natural 
byproduct that may account for the observed structure within our Hubble patch. The data on 
the cosmic background and large scale cosmic structure is rapidly approaching a state where 
this can be tested in detail. I discuss the overall issue, describing and extending post-DMR 
detection work in [l]-[7]. 

1.1. Fluctuation Variables and Their Power Spectra 

Over the length scales we can probe with observations, i.e., within our "Hubble patch" , 
the post-inflation distribution of the noise seems most likely to be linear fluctuations on a slowly 
varying background geometry. 

A generic fluctuation variable D(x,t) can be expanded in terms of modes A4 ^ { adi- 
abatic scalar, isocurvature scalar, vector or tensor; growing or decaying }: 

Di.,t) = /E{«S,w«kA.(x)«kA. + «iS!,*w«kA.w«L} (1) 

f = 1/2 classical , / = 1 quantum . 

For classical fluctuations, aj^-^ is a random variable and aj^.^ its complex conjugate, while 
for quantum fluctuations, aj^-^ is an annihilation operator for the mode kA^ and aj^.^ is 
the creation operator. The mode functions which describe the evolution. The 

spatial dependence of the modes is given by eigenfunctions Qi^j^(^) of the Laplacian of the 
background geometry. For a flat background of most relevance to inflation models, it is simply 
a plane wave, Ok>iW = e'^ ''' labelled by a comoving wavevector k. For curved backgrounds, 
the eigenfunctions are more complex. 

The power spectrum of D associated with mode A4 is the fluctuation variance per log 

wavenumber and can be expressed in terms of the statistics of aj^-^ and aj^.^ : 

quantum: P^l^(fc) ^ ^ = (,)|2 ^ 2(a[^„j^^^) ^ (2) 

classical : ^^[^^(fc) = ^l"k>t(OP (ajkAi^kAi) • (3) 



If the modes are Gaussian-distributed, statistically homogeneous and isotropic, then this is all 
that is needed to specify the patterns in the held D(x,t). 

In the inflation picture, the wavenumbers in the observable regime are usually con- 
sidered to be so high that any pre-inflation mode occupation, (fljj.^flkAi)' negligible, and 
only the unity zero point oscillation term appears. In that case, we connect to the random 
held description by making the real and imaginary parts of aj^-^ Gaussian-distributed with 
variance 1/2. Although quantization is at least self consistent in linear perturbation theory 
about a classical background, there are still obvious subtleties associated with the transition 
from a quantum to a classical random held description. A true inconsistency appears if we 
include the nonlinear backreaction of the fluctuations upon the background flelds and upon 
themselves. For this, we would need a quantum gravity theory. The stochastic inflation theory 
is an attempt to bypass this, by treating the fluctuations quantum-mechanically and the inho- 
mogeneous background classically, with the fluctuations influencing the background through 
stochastic noise terms in a network of Langevin equations for the held variables, e.g., [8, 9]. 
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Over the observable fc-range, it is convenient to separate the issues of overall ampli- 
tude for VD\M{k) — characterized say by VoiMiJ^n^) some normalization length scale R„, 
or, better, by an integral wri a filter, J W(kR„)VD\M(f^)d^T^k, from shape — characterized 
by an index 

noik) + 3 = dlnVD\M{k)/dlnk . (4) 

Thus —no is a "fractal dimension": zero is white noise, while three is scale invariance in D, 
or fiicker noise, with each octave contributing the same loudness. 
In the infiation regime, 

D e {6(t)i„f,6(t)is,h+,h^,6liia,6liiH, 6q, . . .} . (5) 

That is, D would refer to fiuctuations in (1) the infiaton field 64'inj whose equation of state 
can give the negative pressure needed to drive the acceleration, (2) other scalar field degrees 
of freedom b(f)is which can, for example, induce scalar isocurvature perturbations,^ (3) grav- 
itational wave modes h^,h^, (4) the inhomogeneous scale factor a(x,t), Hubble parameter 
iJ (x, t) and 

deceleration parameter : (/(x,^) = —din Ha/ din a , (6) 

encoding scalar metric perturbations and their variations. Infiation ends when q passes from 
negative to positive. Provided the fiuctuations over the observable fc-range remain Gaussian, 
the outcome of infiation is therefore a set of amplitudes for scalar metric (adiabatic) pertur- 
bations, gravity wave modes and various possible isocurvature modes, and primordial spectral 
index functions for each, in particular: 

scalar : ris (A;) = 1 -| ina|ff( ) ^ where (5 In ajij = (5 In a(x, t(x, iJ ~"'^)) , (7) 

(i m /t 

tensor :nt(fc) = -3 + ^^^^^^^^ , where VGwik) = Vh^k) +VhJk) . (8) 

Measuring the power in scalar metric fiuctuations on the time surfaces upon which the inho- 
mogeneous Hubble parameter H(x,t) — the proper time derivative of lna(x,t) — is constant 
is useful [11, 12, 13, 9]: Once Ha exceeds k for a mode with wavenumber k, it becomes time- 
independent during an infiation epoch with a single dynamically-important scalar field, and it 
remains so through reheating and the passage from radiation into matter dominance until Ha 
falls below k (the wave "re-enters" the horizon). ^ 
In the post-infiation period, 

D e {(>Pcdm,^Vcdm,^PB,^VB,(>f~,,(>ferv,^fmv,h+,hx,(^, . . .} . (9) 

That is, D would refer to fiuctuations in the density and velocity of dark matter and baryons 
i^Pcdm , ^Vcdm, 6pB, 6vb), in the distribution functions for photons (Sf-y) and relativistic or 
semi-relativistic neutrinos {6ferv,^fmv), and in the metric (dispersing gravitational wave 

2 If axions are the dark matter, would be the axion field. The isocurvature baryon mode would need to 
have a coupled some way to the baryon number, e.g., [16]. 

3 To be precise about the scalar perturbation quantities used in practice, in the notation of Bardeen [10], 
Slna = 8\nH = H~^^u -\- <I>jj in the 'longitudinal gauge', where H is defined as — l/3x 'trace 
of the extrinsic curvature'. A translation to the time surface on which 8H = gives 51na|jj, which 
[9] used to characterize the metric amplitudes in stochastic inflation. But in linear perturbation theory, 
51na|jj = 5 In a — -^y^-j^SIhH is just Bardeen's ipcom (where — ^com/(4a^) is the 3-curvature on comov- 
ing hypersurfaces) . It is related to Cbst of [H] Cbst — '-Pcom + V2$^/(3(_ffa)2(l + q)). The latter term 
is small when k < Ha: both are nearly constant 'outside the horizon' as long as the inflation models are 
not too outrageous (see § 2.). The use of <-Pcom was advocated by [13, 15] and Cbst by [11, 12, 14]. Either 
will do. 
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modes /i+_x and the gravitational potential for scalar fluctuations $ = —Slna). The Gaus- 
sian nature of the statistics is not modified until mode-mode coupling occurs in the nonlinear 
regime. 

The goal of much of cosmology is to use observations of structure in our Hubble patch 
to piece together the power spectra for observables, then from these to infer the power spectra 
for the post-infiation fiuctuations, i.e., 

{Vir.aURn'),'PGwiR-'),Vs4,jR-'),n,{k),ntik),ni,{k)} , (10) 

and thereby learn about the physics of the early universe. Hampering this program is the 
large number of unknown cosmological parameters. We know well the CMB temperature 
Tcrnh = 2.726 ± 0.005 [17] and the number of light relic neutrinos, hence Sl-^ and, apparently, 
ilerv We do not know the 'global' parameters 

{h, : 7 ^cdm : ^hdm = ^mi^ : ^curv = ^ ^tot : • • •} : (-^-^) 

as well as energy densities and lifetimes for any decaying particles that were once present. (Here 
h is the Hubble parameter in units of 100 km s~^ Mpc~^ and the density parameters are Qj = 
Pj I Per, where p^r = 10.5 h^ kev cm~^.) A reasonably strong case can be made that we actually 
know Sl^h^ = 0.0125 to within 10% or so [18]. The small curvature fiuctuations observed 
with COBE is suggestive of small mean curvature, Q.curv ^ 1- f^A ~ is preferred over the 
odd physics that would be required to make SIa (or (T/((/)))/(3 x 10~"'^^GeV)'*) significant just 
at the current time. The favoured theoretical hypothesis is then that the total density in 
non-relattvtsttc matter, is 1, but, with the best astronomical values for h (~ 0.7 — 0.8), 

one gets a globular cluster age crisis unless SIa or Q.curv is nonzero — or something else exists 
whose energy density varies more slowly than the of nonrelativistic matter. 

As well as the unknowns in 'global' parameters, astrophysical functions required for 
mapping from observable to infiation spectra are unknown. Examples are biasing factors 
relating power spectra for galaxies (Vgg), clusters (Vcc), etc. to those for the underlying mass 
density field (Vpp), 

h]{k) = Vgg(k)/Vpp(k) , hl(k) = VrAh)lVpp{k) . (12) 

The hope is that linear amplification holds over large scales, «,e,,that hg and are k- 
independent, and the power spectra inferred from redshift surveys reveal an underlying density 
spectrum [19]. A prediction is that the power in the cross-correlation of clusters and galaxies 
obeys Vgc(k) = hgh^V pp{k) [20]. Remarkably, the data are roughly consistent with this simple 
picture. 

Another important unknown is the reheating history of the Universe, which may have a 
strong impact upon CMB anisotropies, and because it depends upon when and how efficiently 
massive stars formed in the pregalactic Universe, it is especially hard to predict in a given 
theory. 

1.2. The Observable Range m k-space 

For hierarchical theories of cosmic structure formation, we may roughly divide fc-space 
into various wavebands shown in Fig.l.. (I normalize a to be unity now so that comoving 
wavelengths, 2Trk~^, are expressed in current cosmic length units. Since these are estimated 
from recession velocities, the unit is the h~^Mpc. — 1 is the redshift at time t.) The 

astronomy associated with each band is: ULSS (ultra-large-scale-structure), with k~^ in 
excess of a few times the Hubble radius, cHq^ = 3000 h~^Mpc. We get mean Hubble patch 
values out of this (i.e., "global" parameters such as ^curv) and a little very long wavelength 
fiuctuation information. VLSS (very-large-scale-structure), from the horizon scale (k~^ ~ 
2cHq^ for Einstein-deSitter models) down to say k~^ ~ 100h~^Mpc: 8p and 8v are apparently 
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small enough that they exert little influence on observed cosmic structures, but gravitational 
potential perturbations generate large angle CMB anisotropies. LSS (large scale structure), 
from ~ 100h~^Mpc down to about 5h~^Mpc: we infer that the evolution of the waves in 
this band is sufficiently linear that flrst order perturbation calculations of the large scale 
streaming of galaxies and the clustering of galaxies and clusters may be valid, incomparably 
simpler than trying to correct for complex nonlinearities associated with dynamics and biasing. 
Fig.l. shows the wavebands probed by various large scale structure observations (large scale 
streaming velocities LSSV [21], the angular correlation of galaxies Wgg{9), the power spectrum 
and redshift space correlation function of galaxies as probed by the QDOT and other redshift 
surveys e.g., [22, 23, 25], the correlation function of clusters of galaxies ^^c e,^,, [26, 27, 28]). 
The best indicator for large scale power is the angular correlation function of galaxies [29, 30]. 

CMB anisotropy experiments can be well characterized by fllters which act upon a 
fc-space 'power spectrum for AT/T fluctuations' [31]. Filter functions in fc-space are shown for 
the COBE dmr (~ 7° beam) experiment [32], the firs (3.8° beam) balloon experiment [33], the 
UCSB .sp91 (1.5° beam) 'ACME-HEMT' South Pole experiments [34], the UCSB .sp89 (0.5° 
beam) experiment [35] and the Caltech OVRO ov7 (1.8' beam) experiment [36]. The BigPlate 
Saskatchewan experiment [37] has a similar fllter to sp91, the balloon-borne MAX [38, 39] and 
MSAM [40] experiments have fllters which cover about the same range as sp89. WhiteDish 
[41] and a new OVRO (T beam) experiment cover the region between sp89 and ov7. Of 
more direct observational relevance are the corresponding fllters in f-space shown in Fig.2.(b), 
showing experimental sensitivity to multipole components in the radiation anisotropy pattern. 

Thus, CMB anisotropy experiments cover the entire VLSS and LSS bands. Primary 
amsotroptes of the CMB are those one calculates from linear perturbation theory and which are 
therefore the most important ones because they are easiest to interpret. Their power spectra 
are quite complex [31, 3], because they include effects associated with geometrical ripples in the 
past light cone (Sachs- Wolfe effect), with the flow of electrons at photon decoupling, the degree 
of photon compression at decoupling, and the damping associated with the width of decoupling: 
below ~ 5 h~^Mpc, the primary power is basically erased if hydrogen recombination is standard 
(SR line in Fig.l.); if there is an early injection of energy which ionizes the medium, photon 
decoupling would not have occurred until a lower (Sl^-dependent) redshift and would erase 
AT/T power on scales typically below the NR (no recombination) line shown. 

Below LSS lie wavebands for which gas physics will have been extremely important, 
if not dominant, in determining the nature of the objects we see and how they are clustered. 
Fluctuations are nonlinear in these regimes. The light long-dashed fllter curves at smaller scales 
show the bands probed by very small angle microwave background experiments, the VLA, the 
SCUBA array on the sub-mm telescope JCMT, and the OVRO mm-array. Although their 
beams are too small to see primary CMB anisotropies, they will provide invaluable probes of 
secondary amsotroptes (those generated by nonlinear effects, including redshifted dust emission 
from galaxies and Thomson scattering from nonlinear structures in the pregalactic medium). 

In a hierarchical model, nonlinearity at different scales will occur at sufflciently differ- 
ent epochs that I divide the "gastrophysical" realms into medium, small, very small and ultra 
small, bands ((MSS, SSS, VSSS, USSS), responsible for the construction of, respectively: 
clusters and groups (~ 1{]^'^~^''M.q)] bright galaxies (~ 10"'^"'^""'^^ Mq); dwarf galaxies and Ly- 
man alpha clouds (~ 10®""'^'^ Mq); and the flrst gas clouds to collapse (~ 1O®~^M0), which 
make the flrst stars. Of course, signiflcant gas dynamical processing may obscure the hier- 
archical relationship between object and primordial fluctuation waveband. Further, damping 
processes or tilted initial spectra may require some of the shorter distance structure to arise 
from fragmentation and other non-gravitational effects. 

'Observed' power spectra (actually their square roots) are shown as hatched regions 
for density fluctuations inferred from COBE and for galaxy fluctuations inferred from the 
APM [29] and ROE [30] Wgg data. The long wavelength hatched curve is the rfmr-normalized 
scale invariant spectrum (assuming an = 1 model, and including the current 10% dmr 
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error on overall amplitude). The heavy curve extending the hatched Wgg power into smaller 
distances is the power corresponding to the well known Cgg(r) = (r /rogg)~'' 3D correlation 
function form, where the old CfAl redshift survey values have been taken, r-Qgg = 5.4h~^Mpc 
and 7 = 1.8. Power spectra derived from the QDOT [22], IRAS 1.2 Jansky [23] and CfA2 
[25] redshift surveys are compatible with the range inferred from Wgg when account is taken of 
redshift space distortions and biasing offsets between IRAS and optically identified galaxies. 
As already noted, cluster-cluster correlations and galaxy-cluster cross correlations [28] also 
seem to be compatible with this inferred spectrum. 

The (linear) density fiuctuation power spectra shown in Fig.l. are for three = 1) 

models normalized to the COBE dmr data : a standard CDM model with an initially scale- 
invariant spectra = 1, one with the spectrum tilted to = 0.6, and an = 1 model 
whose shape is characterized by a parameter F = 0.25, whereas F = 0.5 for the standard CDM 
model. To fit the galaxy clustering data requires 0.15 ^ F < 0.3 or 0.2 ^ ris ^ 0.6 (see § 5. for 
discussion). The biasing factor hg is relied upon to move the curves up into the allowed Wgg 
band and nonlinearities to bend the shape upward to match the 1.8 law for < 5h~^Mpc. 
Although this FSS 'extra power' problem has been a subject of intense research on variations 
in the scale-invariant minimal-CDM theme for many years (e.g., [20, 12]), we should bear 
in mind the great success inherent in the extrapolation over so many decades from COBE 
normalization to large and small scale structure formation: it seems scale invariance cannot be 
wildly broken and non-minimality cannot be too extreme, even if the generation mechanism 
has nothing to do with infiation (with the isocurvature baryon model being one deviant case 
[42]). 

2. The Shape and Amphtude of the Primordial Post-Inflation Power Spectra 

During infiation, the same zero point quantum fiuctuation phenomenon which leads 
to the infiaton density perturbations also leads to statistically independent gravitational wave 
perturbations [43, 44]. The equations for the infiaton, isocons, and JJ^^^ ^+,x are derived from 
identical scalar field actions, except the infiaton and isocons are coupled through a potential 
V{4>inj , (t)is, ■ ■ ■), while the gravity waves have no effective mass. Provided the effective masses 
of the scalars are small compared with H^, all respond in basically the same way, rapid oscilla- 
tion of the respective mode functions 'inside the horizon' (k > Ha), almost freeze-out outside 

1/2 

(k < Ha), with a power amplitude {k,t) H/^^tt) essentially given by the Hawking 
temperature on the k = Ha boundary on time surfaces of uniform {Ha)~^, a result moti- 
vated by a WKB treatment of fiuctuation evolution inside the horizon. In stochastic infiation, 
noise at the Hawking temperature radiates from short distances across the decreasing (Ha)~^ 
boundary into a long wavelength background field. We [12, 9] confirmed quantitatively that 
this simple picture, implicit in the early "new infiation" calculations of density perturbations, 
agrees with detailed numerical simulations of fiuctuation generation. The infiaton fiuctuations 
translate into scalar perturbations in the metric, codified in Slna = (H/T1)6(I), where H is the 
background momentum of the infiaton field (f). Thus, the post-infiation spectra are 

pi/^ = y^V^^e-, V!'' =-^^^e-, (13) 
27r ' VqTT mr 2tt ' ^ ' 

The correction factors and Ug to 'the H/(2'it) at k = Ha WKB approximation' are in 
practice nearly zero. How near is now of considerable interest because the COBE results have 
created a desire for calculational precision (§ 2.1.). 

iJ ((/)) and (/((/)) are treated as functions of the infiaton field here, which naturally follows 
from the Hamilton-Jacobi formulation of the problem [9, 45]: the solution to the momentum 

constramt equation, H = fffV is put into the energy constraint equation, turning it 



7 



into the 'reduced Hamilton-Jacobi equation' relating to the potential V((l)): 



with -/?(</)) taking the role of the (reduced) action. (Eq.(14) has corrections dependent upon the 
spatial curvature, hence is valid only for the smoothly varying (long wavelength) background 
field, not the sub-(_ffa)~^ fiuctuating part.) 

The power spectrum ratio and the adiabatic scalar and tensor indices follow: 

nt + 3 r51n_ff-|2_ 1 Vgw 



«t + 3 



1 + g-i +Ct , (16) 



2 

n,,-l 1,-1 -i^^, S^lni? 1,-1 _iSgn{dH/d(f)) m-p dq 

— - — = l + q -q + = 1 + g -q — = j^C,. 

2 At: 0(p^ V 1 + 9 v47r ci<p 

The accurate path to the spectral indices is to take logarithmic derivatives of full numerical 
calculations, a la [12]. The stochastic infiation technique [9] is to write eq.(13) as a function of 
H , q and derivatives, and take a logarithmic derivative wri Ha in place of k, the path adopted 
here and in [9, 2, 3, 5]. Eq.(16) shows that tilt mostly depends upon how far the acceleration 
is below the cn<«ca/ value of unity (but for g — 1, a substantial scalar tilt can come from the 
second term, yet no tensor tilt, as in § 2.4.). Here C't^s are correction factors associated with 
derivatives of the Ui s, which I now discuss in § 2.1., but which the reader may wish to skip 
since I find them to be small and thus drop them subsequently. 

2.1. Corrections to the Stochastic Inflation Calculation of Power Spectra and Their Shapes 
We and others have often used the Ui s = 0,C(_s = approximation e.g., , [2, 3, 5], 
but [44, 15, 47, 48, 46, 49] have stressed the importance of higher order corrections. For the 
case of uniform acceleration, the tensor and scalar equations can be solved analytically in 
terms of Hankel functions and the asymptotic limit can be taken to determine the correction 
factors for the tensor [44] and scalar [15] modes: 



lnr(i-g-i)-g-Mn(-2g)-iln7r , (17) 

-^u[(l)(l + q) + 0(l + qf , 

^-1(1+^1/2...) 
24 



«;(oo) - ^.(1 - ln2) + i ^ \' for large = - q-') 



<{-'!-') ^ ^p^ = -(7 + ln2-l)+[*(i-g-i)-*(^)-ln(-g-i)] (18) 
- -(l-ln2)+ '-"^(l+;;^ + •••) , large.. 

Here ^{v) = d\nV / dv is the diGamma function, 2 — 7 — 2 In 2 = 0.03649 at 3/2, where 7 is 
Euler's constant. The large v limit is surprisingly useful: at —q~^ = 1, it is only off by 2% and 
quickly gets better. There are also weak corrections associated with acceleration changes and 
the effective masses of the scalars, which [46] dealt with by assuming slow changes of these 
quantities to exploit the uniform acceleration analytic solution, not strictly valid but useful to 
indicate the correction level. Following this path, and keeping the leading dq/dcf) term, which 
is the only important one, we have 

u, ^ „, + („;(_5-i) + l)_J__. (19) 
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Thus the correction factors are 

a = 2g-«U-?-^)(l + ?)?-^^^^, (20) 

1 TD^ f)^ In M 

= a+sgn(5F/5,/-)K(-?-') + l)(l + ?)^?-'(5^^^ • (21) 

For evaluations, substituting —0.3 for u'-i.{—q~^) provides enough accuracy (M((1) = —0.27, 
m^(cxd) = —0.31). The key point in Ct is the (l + g) multiplier, which effectively suppresses this 
term relative to the d'^ In H ld(p term of eq.(16): the ratio is 2q~'^u'-i.{—q~^) (l+g) ~ — 0.6(l+(/). 
In § 5.2., we find the data suggests we restrict our attention to tilts < 0.2, hence this ratio 
is below 7%. And when the d^lnH/d^'^ is most important in eq.(16) is when q k, —1, as 
for natural infiation, § 2.4., and in this case, the Ct correction is exponentially suppressed. 
The In 11/8(1)^ correction to Ug has a less strong suppression factor, (1 + q)^^^ , but effective 
enough. An advantage of the forms adopted here over those in [49] is that one is not restricted 
to the (1 + g) regime. But, for the reasons given, I believe it is safe to drop them, which 
I now do. 

2.2. Uniform Acceleration: Exponential Potentials and Extended Inflation 

A constant acceleration regime implies equal scalar and tensor tilts and power law 
infiation (a octP): 

q + 1 = p-\ -l = nt + 3 = -2(p-l)-i . (22) 

Eqs.(14,15) implies an exponential potential, V = Vq exp[— ^167r((; + 1) (ji/m-p]. Of course, q 
must go negative for a viable model of infiation. Nonetheless, over the observable fc-range, the 
exponential approximation is often quite good, even when rather drastic potential surfaces are 
adopted to 'design' spectra. 

Theories with f{<j))R couplings, where R is the curvature, and with one or more dynam- 
ically important scalar fields are a rich source of infiation models. The classical Brans-Dicke 
theory has / = (p /(Abj), where (f) is related to the dilaton. In [12], we considered the induced 
gravity model, with (j) as the infiaton, and showed that if w ~ 10~® in the early universe, the 
coupling of all fields would be effectively weak and the observed density fiuctuation level would 
result. However an arbitrary symmetry breaking potential was invoked to eventually pin 4> 
at (/) = m-p / \/At: to get the observed Newton gravitational 'constant'. In extended infiation 
[50], the infiaton is a separate degree of freedom from the dilaton. The deceleration in the 
conformally transformed frame is uniform, e.g., [51], with value q -\- 1 = 4/(2w + 3), hence 
ris — 1 = + 3 = 8/(2w — 1). Thus another mechanism was required for reheating, bubble 
nucleation. However, to avoid an excessive CMB anisotropy due to large bubbles, the theory 
needed w < 18 — 25 at the end of infiation, yet to satisfy solar system tests to > 500 and hence 
an effective w-pinning or m-p-pinning mechanism is required. 

With conformal transformations, the kinetic term can become nontrivial, making the 
standard Hamilton-Jacobi derivation leading to eq.(15), which predicts q > —1, incorrect. We 
explore the implications of a g < —1 supercritical acceleration for the shape elsewhere [52], 
but for this paper I shall consider only cases for which field reparameterization can take the 
kinetic piece into the standard form, and for which eq.(15) is valid. 

2.3. Slowly Dropping Acceleration: Chaotic Inflation and Power law Potentials 

Power law potentials of the form V^(</') = \em^{(f) / m-p)'^" l{2v) have the advantage 
over exponential laws that q k, — l-|-((/)/m-p)~^i^^/(47r) naturally passes through zero. Chaotic 
infiation discussions [53] have typically focussed on simple potentials, in particular the power 
law form with v taken to be 1 or 2. A characteristic of such potentials is that the range of values 
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of (f) which correspond to all of the large scale structure that we observe is actually remarkably 
small: e.g., for v = 2, the region of the potential curve responsible for the structure between 
the scale of galaxies and the scales up to our current Hubble length is just Am-p ^ </) ^ A.Am-p 
[12]. Consequently, H{<j)) does not evolve by a large factor over the large scale structure region 
and we therefore expect approximate uniformity of ng{k) and nt{k) over the narrow observable 
bands of fc-space, and near-scale-invariance for both. Although this is usually quoted in the 
form of a logarithmic correction to the In ajij-spectrum, a power law approximation is quite 
accurate [2]: 

Ni(k) is the number of e-foldings from the point at which wavenumber k 'crosses the horizon' 
(when k = Ha) and the end of inflation. For waves the size of our current Hubble length we 
have the familiar Ni(k) ~ 60, hence 0.95, 0.97 for i> = 2 and 0.97, 0.98 

for i> = I (massive scalar field case). Further, the observable scales are sufficiently far from 
the reheating scale that Ni is relatively large over the observable range: e, 3,, over the range 
from our Hubble radius down to the galaxy scale, decreases by only about 0.01. 



2.4. Dropping from Nearly-Cnttcal Acceleration: Natural Inflation 

In natural infiation [54, 2], the infiaton for the region of fc-space that we can observe 
is identified with a pseudo-Goldstone boson with a potential V = 2A'* sin^((/)/(2/)). This is 
similar to the axion, except that the symmetry breaking scale / is taken to be of order m-p 
and the energy scale for the potential is taken to be of order the grand unified scale, mauT , 
so that an effective weak coupling, \g = A'^ / ( fm-p)^ ~ (rriGUT / n^v)^ arises 'naturally', giving 
the required 10~^^ for mcjuT = lO^^GeV. To obtain sufficient infiation and a high enough 
post-infiation reheat temperature for baryogenesis, / > 0.3m-p is required. 

To have a tilted spectrum and also get enough infiation in our Hubble patch, 4>/ f 
must have started near the maximum at tt, an infiection point where q is nearly —1 [2]: 

g + 1 ^ (l + ^)exp[-^#z(fc)]R.O, (24) 

r^sik) « + + , nt + 3«0(g + l) , + . (25) 

(."""J ) ' In a|ff 

Thus, we can have a scalar tilt but tensor tilt and gravity wave power are both exponentially- 
suppressed. 



2.5. Rapid Acceleration Changes: Radically- Broken Scale Invariance 

The index can have complex fc-dependent structure when the acceleration changes 
considerably over the fc-band in question. According to eqs.(16), the post-infiation gravita- 
tional wave spectrum will have power increasing with wavelength (the correction C't seems 
unlikely to modify this, although supercritical acceleration can), whereas artfully using the 
d^/l + q/dcf) term in the infiaton effective potential allows essentially any prescribed shape for 
the adiabatic scalar spectrum. 

A prion, it seems unlikely that a marked change in the expansion rate or acceleration 
would just happen to be in the narrow window of fc-space accessible to our observations. How- 
ever, in (/)-space, this window is not very far from (fignd defining the acceleration/deceleration 
boundary, hence the q rise to zero must be reasonably rapid in (f). Even so, for the models 
described above, the rapid change does indeed occur only near the end, suggesting special 
physics might have to be built in. 
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Rapid acceleration changes, if present, would seem to be more likely a consequence of 
interaction with other field degrees of freedom rather than a result of infiaton self-interaction. 
Thus, many of the toy models constructed to illustrate that radically broken scale invariance is 
possible involved two scalars interacting with either simple polynomial potentials (with second, 
third and fourth order terms) [55, 12], or combinations of exponential potentials [47]. 

Even with many scalar fields being dynamically important, it is often possible to con- 
sider an effective single infiaton self-interacting through an effective single-infiaton potential 
over the observable scales. This is because the fields first settle into gorges on the potential 
surface, then follow the gorge downward towards the local minimum along a single field de- 
gree of freedom, 4>\\, to be identified with the infiaton. The other degrees of freedom, 4>_\_, are 
'isocurvature' degrees of freedom. Usually, the faces rising up from the gorge will be suffi- 
ciently steep that the inevitable quantum noise that excites motion up the walls quickly falls 
back, leaving no usable isocurvature imprint, effectively making those dimensions irrelevant 
(although curvature in the trough can lead to complications in the kinetic energy piece of the 
infiaton degree of freedom). 

Many models of double infiation could be described this way, consisting of two periods 
of infiation with relatively constant H , one at high H , the other at low H . These lead to nearly 
scale invariant fiuctuations in the two associated regions of fc-space, high amplitude, then low. 
The join must be accompanied by a large change in acceleration, hence in Ug ti^) over the 
corresponding fc-band: exactly how one crafts the transition determines the detailed form of 
rig tik). General variations of the effective single infiaton potential H{<j)), hence of V{(j)), allow 
wide latitude in what can be constructed. Since Ug has a term cx -\- q/d(f), it tends to be 
more susceptible to the variations than rit is, and therefore the adiabatic scalar spectra should 
exhibit sharper structural features. 

Models with two scalar fields that do not allow an effective single infiaton approx- 
imation over the relevant band in fc-space have also been used to construct power spectra 
with mountains and valleys and also to generate non-Gaussian infiation fiuctuations. Often 
these involve an instability, with negative transverse components of the mass-squared matrix, 
d'^V / d(f)id(f)j , leading to an opening up of the gorge or its bifurcation. Tuning the location of 
such a structure to the window on the potential surface we can access must be unpalatably 
precise [12]. 



3. Inflation-Based CMB Power Spectra c.f. the Data 



3.1. Theoretical CMB Power Spectra 

For a given infiation model, perturbed Einstein-Boltzmann equations {e.g., [56, 31, 3]) 
must be solved for each mode Ai present to get the temperature radiation pattern at our 
location and at this time: 

[q, here, now) = ^ = }^ 'Yim{q) , (26) 

where /y is the unperturbed Planck distribution and Sf^'^^ is the distribution function fiuc- 
tuation. If the post-infiation fiuctuations are Gaussian-distributed, then so are the multipole 
coefficients a^^"*, with amplitudes fully determined by just the angular power spectra C^^"*, 
defined by 

(27) 

Ztt 

Sample theoretical C^'s are shown in Fig. 2. (a). The "standard" adiabatic CDM model 
(SI = 1, ris = 1, h = 0.5, Sis = 0.05) with normal recombination illustrates the typical 
form: the Sachs- Wolfe effect dominating at low £, followed by rises and falls in the first and 
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subsequent Doppler peaks, with an overall decline due to destructive interference across the 
photon decoupling surface. A similar CDM model, but with early reionization (at z > 200), 
shows the Doppler peaks are damped, a result of destructive interference from forward and 
backward flows across the decoupling region, illustrating that the "short-wavelength" part 
of the density power spectrum can have a dramatic effect upon Ct, since it determines how 
copious UV production from early stars was. Lower redshifts of reionization still maintain a 
Doppler peak, but suppressed relative to the standard CDM case. 

A form often adopted to describe the low-f end is the "Sachs- Wolfe" power for scalar 
metric perturbations [31] 

= I f d\nkV^{k)jj{kRr) , (28) 

^ rff-M A^-- r(£-i^)r(£ + 2) r(i-i^)r(|) 

~ g/^i^n n^j r(£)r(£ + 2 + i^) r(f + i^) ' ^ ' 

(l-n,)(ll + 6(l-n,) + (l-"0') 
' - 24 • 

Here is the comoving distance from us to the surface over which photons decouple from 
the baryons, about 5800h~^Mpc away for a CDM cosmology, and i?„ is a normalization scale. 
The analytic result in terms of Gamma functions holds if there are no deviations from the 
power law form for the power spectrum of the gravitational potential V<^, related to V\na\H 
by ~ (5/3)'P^^^. This predicts a flat Ci for rig = 1, obeying the pleasing formula 

C]''^ = K ■Pil/fi^/S, relating AT/T to $/3. In fact, for a realistic model, there are 

corrections to this formula from other anisotropy sources: in particular for the standard CDM 
model shown, there is a small rise over the multipoles that COBE probes, modelled by an 
effective index 1.15 if we use the eq.(29) form. 

The dotted Ci in Fig. 2. (a) also has a flat initial spectrum, but has a large nonzero 
cosmological constant in order to have a high Hq, in better accord with most observational 
determinations. The speciflc model has SIa = 0.75, Slcdm = 0.22, Sl^ = 0.03, -ffo = 75, = 1. 
As one goes from £ = 2 to £ = 3 and above there is flrst a drop in Ci [57], a consequence of the 
time dependence of $ which results in corrections to ~ V^/^ /3. 

A major goal of CMB anisotropy research is to determine all of the ups and downs 
of the spectrum in detail. The 'data points' of Fig. 2. (a) are averages of the C^'s wri 'fllter' 
functions Wi'. 

{Ci)„ = l{CtWt\ll{Wt\ , where = J2 ^^^^ (3°) 

deflnes the "logarithmic integral" 2[fi\ of a function fi. The Wi are taken to be those for a 
set of existing anisotropy experiments spanning the range 10° to 2' shown in Fig.2.(b): (Ci)-^ 
then characterizes the broad-band power that the experiment is sensitive to. The location in 
f-space is {£)-^. The error bars are 10%, the best dmr can possibly do with the full 4 years of 
data, and ones that are actually quite reasonable for intermediate- and small-angle mapping 
experiments. The observed data points for these experiments are shown in Fig. 4. below. 

If the spectrum is tilted, gravitational waves will generally be present to induce a 
tensor-mode spectrum, C^*-* to add to the adiabatic scalar spectrum C'f \ The amplitude of 
gravitational wave modes decays by directional dispersion as the modes re-enter the horizon, 
just as waves in any relativistic collisionless matter do. Before the gravitational waves disperse 
however, they influence the microwave background through the inhomogeneous redshifts they 
induce. The recognition of the potential importance of gravity waves for large angle microwave 
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background fluctuations has a long history, e.g., [43, 44] and generated a tremendous post- 
COBE burst of excitement, and papers e.g., [61, 47, 48, 62, 2], especially when the sensitivity 
to tilt was realized that meant the rfmr-signal could even be largely tensor-induced. I shall 
parameterize the relative magnitudes of scalar and tensor by the ratio of the broad-band powers 
associated with the rfmr-beam's fllter: 

\ t I dmr 

This ratio has no simple analytic result and its value is dependent upon the details of the 
cosmology being considered. The 1.5 numerical result holds for small deviations from scale 
invariance, Ug k, rit, and for CDM-like models. (In [3, 5], we used instead C2*V^2^'' ~ l-2?ts to 
characterize the relative magnitudes.) 

In Fig. 2. (a), the dashed line shows a standard CDM model, but with a chaotic-inflation 
inspired rig = 0.95 tilt, along with a fts ~ 0.3 gravity wave contribution (eq.(31) with rit = rig), 
contributing the lower dashed line to the total [3]. Although d"p has a flat part at low I just 
as df ^ does for this nearly scale invariant model, there is about a 20% drop from £ = 2 to 
£ = 3, and there is no Doppler peak, only a rapid decline at £ > 50. 

To get a visual impression of what the spectral structure means. Fig. 3. shows what the 
sky looks like on a few resolution scales for the standard rig = I CDM model: on the COBE 
beamscale (Gaussian flltering £g = 19), the nearly scale invariant form; on the half-degree scale 
(£g = 269), where the standard recombination spectrum is a maximum; with no smoothing at 
all, with the shapes deflned entirely by the destructive interference that occurred across the 
photon decoupling region. For early-reionization, the shapes in the 60° NR map are also the 
naturally occurring ones, since there is no power left at £g ~ 269 to artiflcially fllter. 



3.2. How Accurately Can The Spectra Be Measured? 

In [5], we showed that for small variations about the CDM model, the height of the 
flrst Doppler peak relative to the dmr band-power is (within ~ 15%) 

Ci\^ ^ ^ ^-3.6(i-n,) ^ rig-Hl + ru)/3.6-0.b[(l-n^)ih-^^]-(^^y" , (32) 

where Qb = 0.0125h~^ has been flxed at the standard BBN value. Here zr is the reionization 
redshift and must be < 150 to have a local maximum. For example, a model with no gravity 
wave contribution (as natural inflation would predict) but rig 0.88 has a spectrum that 
is almost degenerate with the rig = 0.95, fjs = 0.3 spectrum, so much so that it will be 
extremely difficult to tell them apart. More generally, we argued that the precision required to 
separately determine rig, rig, SIa, ... is too high for likely experiments, but rig can determined 
accurately. (An exception is Sl^h^, which the relative heights of the Doppler peaks are sensitive 
to.) To separate the various contributions to rig requires other cosmological experiments, 
e, 3,, measuring the scalar perturbation shape through galaxy-galaxy power spectra (§ 5.1.) 
and amplitude through cluster abundances or streaming velocities (§ 5.2.); and, in some happy 
future, determining Hq definitively. 

We now discuss the experiments in more detail, flrst in an idealized way to show what 
is needed for achieving even the 10% error bars shown on the band-powers in Fig. 2. (a). The 
signal (AT/T)p from the p*'' pixel of a CMB anisotropy experiment can be expressed in terms 
of linear fllters J^p/m acting on the a^™: (AT/T)^ = Y^i^^p,imaim- The Tp^im encodes 
the experimental beam and the switching strategy that deflnes the temperature difference, the 
former flltering high £, the latter low £. The pixel-pixel correlation function of the temperature 
differences can be expressed in terms of a quadratic Npi^ x Npi^ fllter matrix Wppi / acting on 
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Cti'- 

'^Tpp' = ((~) (~) ) = l[Wpp'/CTf^ , Wppi^t = 2£_^ ^ ^^P,lm^p',tm (33) 

Im 

= (— ) =-—Y,Ctpv=^ \WtCTt\ , "^^^ = E ^PP.^ • (34) 



p = l '- p = l 



The trace Wi defines the average filters [58, 59, 60] shown in Fig.2.(b), which determine 
the rms anisotropies (Tt[1/F]. Typically we will be given the anisotropy data in the form 
(AT/T)p ±aop, where aop is the variance about the mean for the measurements. In general, 
there may be pixel-pixel correlations in the noise, defining a correlation matrix C'opp' with 
off-diagonal components as well as the diagonal (r'jjp. 

In the simplest experiment that can be imagined, we would have Copp' = Cd^pp' ^^'^ 
the pixels sufficiently separated on the sky that only Wi is an effective probe of Cf, «,e,,that 
Ctpp' ~ ^T^pp'- {^i) B th = '^T I'^i^ i) is the quantity we wish to estimate. For large Npi^, 
the 1-sigma uncertainty in the experimental value of the band-power is [5, 6, 7] 

{Cl)B.ois = {Cl)B.ma.L±\/m7^[{Cl)B.ma.L+^lmW^]]; (35) 
{Ci)B,ma.L = {Ci)B,ik ± ^/W^^ [{Ci)B,ik + ^^1/^1]] ■ (36) 

An experimental noise uu below 10~® is standard now, and a few times 10~® is soon achievable; 
hence, if systematic errors and unwanted signals can be eliminated, the 1-sigma relative uncer- 
tainty in {Ci) g will be from cosmic-variance alone, ^J^/Nd, falling below 10% for Npi^ = 200, 
I.e., a mapping experiment. For large Npi^, the observed maximum likelihood will fiuctuate 
from {Ci)g the quantity we want, according to eq.(36), but the error bars of eq.(35) include 
these realization-to-realization fiuctuations (thus \/2 appears, not 1). If there were full-sky 
coverage and errors from cosmic variance alone, the fractional error in the {Ci)^, goes as 
~ (£) ^ [7]: so tiny for intermediate and small angle experiments that it would appear that 
even extremely subtle differences in the spectra could in principle be determined at high I. 
Now I shall discuss how we are doing so far in practice. 

3.2.1. Current Status: Experimental Broad-Band Powers 

To determine band-powers for an experiment [4, 7], I construct a local model of Ci, 
assumed to be valid over the scale of the experiment's average filter W f. 
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{Ci)b (^ + i)'+""" I[Wi]II[Wi{l + i)2+"--] , « - 3 . (37) 



The local "spectral colour index", mat, is similar to the nu of eq.(4), except for 2D, not 3D. 
The form differs very little from that in eq.(29): for small I, n^T is related to rig as shown, but 
with rig now interpreted as a phenomenological rather than a primordial index. Thus = 3 
corresponds to white noise in AT. To get the band-power, I use eq.(37) with n^T = —2, 
i.e., scale-invariant over the band, but check that the result is robust to variations in mat- 
This is true for all intermediate and small angle experiments to date, and as we shall see, even 
holds true for dmr and firs, which have such a large coverage in f-space that they can also 

1/2 

be used to determine the index mat- The amplitude (Ci) g can be determined by whatever 
statistical method we are most enamoured with, whether Bayesian as I prefer, or frequentist. 

There are so many detections now that I split Fig. 4. into two panels for clarity, the 
upper giving the overview, the lower focusing on the crucial first Doppler peak region. These 
figures have been evolving rapidly since I introduced them [4, 7]. Data points either denote 
the maximum likelihood values for the band-power and the error bars give the 16% and 84% 
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Bayesian probability values (corresponding to ±1(7 if the probability distributions were Gaus- 
sian) or are my translations of the averages and errors given by the experimental groups to this 
language. Upper and lower triangles denote 95% confidence limits unless otherwise stated. See 
[7] for details. The horizontal location is at {£)^r and the horizontal error bars denote where the 
filters have fallen to e"*^^ of the maximum (with Fig.2.(b) providing a more complete represen- 
tation of sensitivity as a function off). Only wavelength-independent Gaussian anisotropies in 
AT/T are assumed to be contributing to the signals, but non-primary sources (e.g., dust, syn- 
chrotron) may contribute to these C^'s (as can unknown systematic errors of course). Cleaning 
AT/T observations has been done to some extent in most of these cases, and will be the key to 
the ultimate accuracy that we can achieve in spectrum determination - not theoretical cosmic 
variance. Generally the underlying primary Ci will be lower than the values shown, but it can 
be higher because of 'destructive interference' among component signals. 

Proceeding from small £, the £ = 2 power uses the first year 53GHz quadrupole values 
with a Galactic cut hocut = 20° [4]. The value for the combined first and second year data 
will be even lower [63], but it is also the multipole most likely to have a residual Galactic 
signal contaminating it, possibly destructively. The solid dmr point with the tiny error bar 
is my translation of the combined first and second year result [63], while the open dmr point 
for the first year data and the firs point are my band-powers, from refs.[4, 7] and § 4.. The 
Tenerife point uses the data for the limited region of the sky they probed at both 15 and 33 
GHz. Remarkably, in view of the relatively low frequency, the band-power for their data at 
15 GHz only, which covered a much larger region of the sky, agrees. We now come to the 
confused region from two degrees to half a degree, which can be better seen in the lower panel. 
Two sp91 band-powers are shown, for a 9 point scan and a 13 point scan [34]. All 4 channels 
were simultaneously analyzed [4]. The offsets are for clarity. The BigPlate result [37], bp, 
is slightly higher. New sp94 results apparently also have detections at a little higher level 
and there will also shortly be a new bp result, both with more extended frequency coverage. 
Python [65], py, has wide coverage in f-space, but has only the single 90 GHz frequency so 
signal cleaning cannot be done. Argo [66], ar, is next. The next 3 results are for 3 scans from 
the fourth fiight of the MAX balloon-borne experiment [39], M4, the open squares for the Iota 
Draconis scan and the Sigma Hercules scan, the solid point for a GUM scan. Because the 
filters changed with frequency, the points are placed at the average over all MAX4 filters. The 
two solid data points are for the third MAX fiight [38], the upper for a GUM scan, the lower 
for a Mu Pegasus scan (which had a strong dust signal removed). The GUM point also shows 
the Bayesian 1 — sigma error bar. The dotted lines ending in triangles denote the 90% limits 
for the MSAM [40] single (g2) and double (gS) difference configurations for all of their data, 
although there are some worries that half of their data was contaminated by non-Gaussian 
sources, which, when excluded, lowers the band-powers somewhat. The next three give upper 
limits, but no detection. The open triangle is the 95% credible limit for the sp89 9 point scan 
[35, 60]. Switching back to the upper panel, the 95% limit from the m = 2 mode analysis 
of the WhiteDish experiment [41] is wd2. There is a hint of a detection in the m = 1 mode 
analysis. The ovro 7 point upper limit [36] is last. New ovro experiments with significantly 
higher sensitivity for this filter and for the ovro22 one in Fig.2.(b) also have results that will 
soon be available, once the detections have been cleaned of radio sources. 

In the future we will be able to strongly select the preferred theories by simultaneously 
analyzing experiments like these. Although we have already tried this, e.g.,ioi sp89 and 
ov7 [60] and for dmr, sp91, sp89 and ov7 [3], for the time being I believe it can be quite 
misleading because we cannot be confident yet that the data has been properly cleaned of 
secondary backgrounds, foregrounds and instrumental systematics to reveal the underlying 
primary anisotropies. Until then, band-power figures such as Fig. 4. should be our guide to 
the evolving progress towards a primary Ci spectrum, and the theory of fiuctuation generation 
underlying it. 
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4. DMR and FIRS 

In this section, I describe some Bayesian results on the firs and first year dmr maps 
which use all aspects of the maps simultaneously and so are highly sensitive to all components 
in it. Since we wish to extract only the primary signals, filtering unwanted residual signals is 
essential, but there is some danger in doing so. I compare my unfiltered and filtered results 
with those obtained by the dmr and firs teams, whose techniques also employ filtering of one 
sort or another. 

4.1. Stgnal-to-notse Etgenmodes for Maps 

A full Bayesian analysis of maps requires frequent inversion and determinant evalua- 
tions of Npix X Npix correlation matrices, the sum of all C'Tpp' in the theoretical modelling plus 
the pixel-pixel observational error matrix C'opp'- To facilitate this, I expand the pixel values 
(AT/T)p into a basis of "signal-to-noise" eigenmodes for the maps in which the transformed 
noise and transformed (wanted) theoretical signal we are testing for do not have mode-mode 
correlations, «,e,, which are orthogonal. This can always be done, no matter what the experi- 
ment. Complications are associated with the removal of averages, dipoles etc. and the existence 
of secondary signals in the data, both of which do couple the modes. A model for the various 
contributions that make up the observed data is then 

p=i 

{ukUk') = Skk' , {skSk') = ^TRkhk' 

(resfcresfc') = TZkk' , {ckCk') = ICkk' ■ 

Because the transformation Cj^^^^CtCo^^'^ has dimension (theory variance/ pixel error^), I 
call the £TR,k S'/A^-eigenvalues. I sort the modes in order of decreasing £TR,k, so low fc-modes 
probe the theory in question best. 

Because there are an equal number of eigenmodes as pixels, this new expansion is a 
complete (unfiltered) representation of the map. With uniform weighting and all-sky coverage, 
the eigenmodes are just the independent Re(aim) and Im(aim), with the lowest £ having the 
highest £TR,k, hence k ^ (£ + 1)^. With Galactic cuts followed by dipole removals, and 
especially with inhomogeneous pixel coverage - a bigger issue for firs, but important for dmr 
too - they are complicated and theory-dependent. However, the high S/N-modes are indeed 
the ones that involve large scale pixel linear combinations, while the low S/N-modes typically 
involve positive and negative contributions from nearby pixels that are not sensitive to large 
scale structure in the maps, but are quite sensitive to physics inside the beam, whether from 
systematic effects or true white noise on the sky. This suggests this can be an ideal set for 
filtering. Filtering using S/N-modes has a long history in signal processing where it is called 
the Karhunen-Loeve method [67]. 

There is an arbitrary average and dipole that can be added to the maps, Ck, which I 
take to be described by a Gaussian with very wide width, i.e., a uniform prior probability. In 
S/N-space, this contribution has off-diagonal correlations which affects small k. The residuals 
are modelled by an excess pixel noise with an amplification factor r, which soaks up a significant 
part of the excess I observe in the data, and an unknown component denoted by res^. Of course 
without identifying it, we do not know its correlation matrix TZ, but the data itself can tell us 
something about its structure. 

The sum of \£,k\^ over bands in S/N-space defines a S/N power spectrum which is 
easy to interpret because the modes are basically independent of each other, but have the 
disadvantage of depending upon the theory being tested for. In [7], I showed how dmr and 



= Sk r)nk +Ck+ lesk , k =1,. . ., #p438) 

(39) 



Table 1: {Ci)]^^^/!^ ® as a function of rfmr map ( 


"AT + 3 = 


1) 


53a+b 


53a+b(7°) 


53a-b 90a+b 


90a-b 


31a+b 


firs 


1.03±.15 


0.99±.14 


0.30±.30 1.08±.21 


0.00±.30 


0.89±.35 


1.09±.26 




Table 2: {Ci) 


dmr, fir. /lO ^ ^s a function of n, E 


= n^T + 3 




1 n 
o .u 




2.0 1.5 


1.0 


0.5 


0.0 


1 1 0-1- 1 n 


1 n7-i- 1 1 

1 .U ( It . 1 1 


rfmr53a+b 
1.05±.12 1.02±.14 


1.03±.15 


1.05±.19 


1.07±.22 


1 .oozt .zo 


1 .ooit .zo 


firs 

1.24±.25 1.15±.26 


1.09±.26 


1.05±.29 


1.04±.32 




Table 3: 


: (C.>l'r,f,r./10-'as 


a function of kc^t 






142 


132 102 


72 


52 


42 


1.03±.15 


1.02±.15 


rfmr53a+b, n^T + 
1.02±.15 0.94±.14 


■3 = 1 
0.85±.14 


0.87±.18 


0.98±.24 


1.06±.12 


1.05±.12 


rfmr53a+b, n^T + 
1.05±.12 0.97±.12 


■3 = 2 
1.02±.15 


1.05±.21 


1.21±.29 


1.11±.27 


1.07±.26 


firs, riAT + 3 = 
1.05±.26 1.00±.24 


= 1 
0.99±.25 


0.89±.26 


1.01±.32 


1.27±.24 


1.21±.23 


firs, riAT + 3 = 
1.17±.23 1.08±.23 


= 2 
1.12±.27 


1.14±.32 


1.14±.38 


dmrl: S92 
0.97±.28 


Table 4: {Ci)]^^^ firs/^^ ^ from the dmr and ^rs teams 
dmrl: W94 dmr2: B94 dmr2: B94 (mat + 3 = 1.59) 
0.97±.16 l.OOi.lO 1.02lj? 


firs: G94 
1.08±0.3 



firs spectra for both n^T +3 = 1 and 2 reveal excess power at low S/N (high k) in the data, 
which neither beam-filtered theory can account for. Adding the constant r ^ pixel noise 
enhancement is quite a good model for the excess S/N power at the high k end of the firs 
data, and not as good a model for the dmr 53a+b GHz data, but for dmr the map-dependent 
preferred r is significantly smaller than for firs. For both S/N power spectra, there is an excess 
at k about I42 that the power law theories cannot account for, and that plagues the statistical 
analysis of mat- The nature of the modes and of the residual can also be probed by testing 
for its angular structure with correlation functions for S'/A^-filtered maps with high pass and 
low pass filters (as described in Fig. 7.). 

The first step in the Bayesian method is the construction of a joint likelihood function 

i° (^^)dmr,/irs ' "AT and r. For given mat, these reveal a strong maximum in (Ci:>dmr,/irs • 
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Integrating over r (marginalizing it) allows one to construct n^T-{Ci) ^^j^,. j^^.^ contour maps. 
Results for the 50% Bayesian probability value of {Ci)^^^,, with 'one-sigma' error bars de- 
termined by using the 84% and 16% Bayesian values are given in Tables 1-3: Table 1 gives 
band-powers for fixed mat -|-3 = 1 for various map combinations, which show very good agree- 
ment between the 53, 90 and 31 GHz dmr a+h maps and the firs map, and with no discernable 
signal in the dmr a— b maps; Table 2 shows that the derived band-power is very insensitive 
to the index mat! Table 3 shows (for mat -|- 3 = 1 and 2) that the derived band-power is 
remarkably insensitive to the signal-to-noise cut. kc^t gives the number of modes kept, to be 
compared with 928 modes for the dmr maps (5.2° pixel size chosen. Galactic latitude cut of 
25°) and 1070 for the firs map (2.6° pixel size chosen). In Table 4, I compare my results with 
the original first year dmr number [32], a later update [68], the new result using the combined 
first and second year data [63], derived for mat -|- 3 = 1 and also for the most probable index, 
and the recent firs team result using the correlation function [69]. 
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Table 5: 


«AT + 3 for rfmr53a+b and firs 




dmr dmr(7°) 


dmr, A; < (14)^ dmr, k < {U f 


firs 


^•'-'-0.4;-1.0 ' -0.4;-0.9 


1 r+0.6; + 1.4 , „ + l.l; + 1.6 
-'-•'J-0.7;-1.0 -'-•'-'-0.7;-1.0 


-'-•°-0.8;-1.3 



Table 6: mat + 


3 dmr Team, First Year, 


and firs Team 


dmrl: S92(CF) 


dmrl: W94(PS) 


firs(CF) 


"J-0. 65 


1 RO+O 45 


-'-•'-'-0.5 



Table 7: mat + 3 


dmr Team, Secoii 


id Year Indices 




dmr2: B94(CF) dmr2: 


W94(PS,3-19) 


dmr2: W94(PS,3-30) 


dmr2: G94 (3-30fil) 




-'-•^'^-0.44 


-'-•^'J-0.45 





I conclude that the band-power is very robust and well determined, at nearly the same 
level for all of the dmr maps and for the firs map, and largely independent of rig. By contrast, 
the quadrupole power C2 (or equivalently Qrms,PS = 2.726A' ^/ (5/12)^2) varies considerably, 
being related to the band-power by 

^ ^ \ 17 6^K ) ^ ^'^^1 dmr ' ttdmr = 0.73, a/i„ = 1.1(40) 

1/2 

The anisotropy colour index mat = «s— 3 is another story. Integrating over {Ct) j^^^ 
as well as r gives the probability distribution for mat- Table 5 gives colour indices for 
rfmr53a-|-b, both unfiltered and with k^ut in the region where the S/N power spectrum re- 
veals the excess that single mat laws can't explain. Maximum likelihood values with Bayesian 
1-sigma and 2-sigma errors are shown: the index is determined with significantly less precision 
than the band powers are. The k < (16)^ result is similar to the unfiltered result. As expected, 
sharp S'/A^-filtering does lower the index as the residual seen in the S/N power spectra plots 
in [7] is chopped off. The signal-to-noise for modes with k > 13^ is below 0.3, whereas the 
lowest 50 dmr and 20 firs modes have signal-to-noise in excess of unity. The non-diagonal 
part of C'd does not change the Bayesian 50% value in Table 1, and increases the most likely 
band-power by 1.6%. The calculations used the rfmr beam and digitization and pixelization 
corrections advocated in [70]. To check that a faulty beam is not the problem, I tried a smaller 
7° beam with no pixelization corrections, which lowers the value by 0.3. (The corresponding 
band power is listed in Table 1). A physical mat = white noise source as well as a long 
wavelength mat ~ —2 would work quite well, but I haven't explored it yet, although allowing 
r to fioat as a function of cut mimics the effect: that strategy robs the scale-invariant band- 
power by ~ 20% to feed an increased r [7]. Tables 6 and 7 give various results obtained by 
the dmr and firs teams, using correlation function (CF) analysis (S92 [32], B94 [63] and firs 
[69]), quadratic power spectrum estimation (PS, W94 [71], with index determination over the 
£ range shown) and using a linear multipole filtering method (G94 [72]). 

4.2. Sharp S/N -filtering and Wemer-filiermg of Maps 

I have found sharp S'/A^-filtering preferable to smooth S'/A^-filtering for statistical 
analysis, but some examples of the statistical use of smooth S'/A^-filtering were given in [4, 7]. 
One immediate byproduct of the S/N eigenmode expansion is optimal or Weiner filtering, 
very useful for constructing maps cleaned of noise (Figs. 5. ,6.) to show robustness of structure 
from map to map. In Ref. [73], Bunn ei a/, have independently applied Weiner filtering, to 
the 'reduced galaxy' dmr map, a linear combination of the 3 frequency maps. Because the 
Weiner filter changes with both mat and {Ci) ^^^.^ , I have not found it to be useful for the 
statistical analysis of either band-power or spectrum-colour. 
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Given observations ^k, the mean value and variance matrix of the desired signal Sk 
are (see e, 3,, Appendix C of [19]) 

= y^\STR[STR + {l + rf +TZ + IC]~'} 6' (41) 

^ — ' J kk' 

k' 

(As, As,, 10 = \£TR[STR + {l + rf +TZ + IC]~'[{l + rf +TZ + IC]} . (42) 

J kk' 

The mean field (s, is the maximum entropy solution. The operator multiplying is the 

Weiner filter. With no constraints and no extra residual, it is just £TR,k/{{i + r)^ + £TR,k)- 

The fluctuation of the signal about the mean is As, = s, — (s, j^, a realization of which is 

found by multiplying a vector of Npi^ independent Gaussian random numbers by the square 

1/2 

root of the variance matrix {[STR,k/{^ + STR,k)] (1 + r) with no constraints). 

The Weiner filter depends upon the overall amplitude we think the theory signal has, 
parameterized by {Ci) ji^^ ■ I use the maximum likelihood values of {Ci) ji^^ in the 
following. When the noise is large, as it is for these maps, it is the higher I power that is 
preferentially removed by optimal filtering. Thus the theoretical fluctuation As, would have 
to be added to the Weiner-filtered maps of Figs. 5. ,6. to have a realistic picture of the sky given 
the data and the theory. That is the maps are too smooth, and more so for the noisier 90 
GHz map than the 53 GHz map. Still it is very encouraging that the same large scale contour 
features persist in both maps. In the even noisier 31 GHz map, only a hint of the features 
persist visually because the S/N filtering is so strong, but it is evident in the correlation 
function (Fig. 8.). 

The distinguishing feature of the filtering procedures used here and in [72] is that 
they act linearly on the pixel amplitudes. Although using combinations, (AT/T)p(AT/T)pi , 
of pixel pairs is more complicated for statistical analysis, very useful forms of filtering become 
straightforward with these quadratic statistics. The correlation function is the classic example. 
The pixel-error enhancement model for the residual (parameter r) only contributes to the zero 
angle bin of the correlation function. The excess residual (beyond pixel-enhancement) resides 
between k about 16^ to 13^. Fig. 7. (a) shows the dmr correlation function for modes above 
and below these values. Outside of the ~ 6fwhm beamsize, C'(6) is nearly zero for the low S/N 
modes and is nearly the same as C'(6) with no cut for the high S/N modes. The same story 
basically holds true for firs, as can be seen in (d), except the deviation extends to ~ 29j^hm- 
(Note that what is uncorrelated in the full map should have some residual extended correlation 
in the filtered map). The other panels show the correlation function for a map realization of 
pure signal (b,e) and pure noise (c,f). For the pure noise simulation, the cut does matter, 
while for the pure theory simulation, £TRk is so small above 16^ that it doesn't. I conclude 
from this that because C{9) filters small-angle systematic (or physical) effects (if the first few 
angular bins are downweighted), n^T determined this way provides a good indicator of the 
angular colour of the large-angle sky. Although bins within Ojwhm are usually included, they 
are not overly weighted relative to larger angle bins. 

In Fig. (8.), I contrast the highly noisy correlation functions (showing the error bars) 
for the unfiltered maps with C{9) obtained for the Weiner-filtered maps. They reveal a num- 
ber of remarkable points: (1) the optimally filtered correlation functions are effectively mat- 
independent (shown are the n^T + 3 = 1,2 cases for rfmr53a-|-b and firs in (a) and (c)). Of 
course this does not address how statistically significant the optimally-filtered correlation func- 
tion is for a given value of mat- (2) The optimally filtered a-b rfmr maps have zero C{9). (3) 
The optimally-filtered a-|-b dmr C{9ys shown in (f) on an expanded scale agree rather nicely. 
This is especially gratifying for the very noisy 31a-|-b map. The depth of the negative trough 
at ~ 50° damps with increasing map noise, consistent with Weiner-filtering getting rid of the 
higher I contributions. The next step is to apply the S/N eigenmode method to the combined 
first and second year dmr data. This should be better behaved statistically, although [72] 
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point to a continuing problem in the rfmr53a+b map for 14 < £ < 19 modes, just where I have 
found anomalies that cannot be fit by a simple power law that are largely responsible for the 
steep indices. 

5. Large Scale Structure Constraints 

5.1. Parameienzmg the Shape of the Density Power Spectrum 

Shape constraints derived from large scale structure observations rely on the assump- 
tion of a uniform biasing and linearity on large scales, with different biasing parameters for 
the different objects observed. If so, the galaxy-galaxy and cluster-cluster power spectra and 
correlation functions directly reveal the underlying linear density power spectrum, Vp(k). It 
is useful to discuss these observations in terms of physically-motivated parameters that define 
the curvature of Vp(k), the initial scalar spectral index ns(k) and the scale of the horizon at 
redshift Q„r/^er when the density in nonrelativistic matter, Q„r(i~^, equals that in relativistic 
matter, Qer(i~^, 

k'l^ = 5 h-^Mpc , r = fi„rh [fieJ(1.69fi^)]-i/2 ^ (43) 
In [1], we adopted the now- widely-used functional form: 

Vp(k) (X k^+"'(''^ \^l + [ak + (hkf'^ + (ckfYy^'\ (44) 
(a,&,c) = (6.4,3.0,1.7)r-^h-^Mpc,iy= 1.13. 

For r = 0.5, this accurately fits the linear power-spectrum of the standard CDM model with 
= 1, h = 0.5 and Sl^ 0.05 [56]. The oft-used Sl^ form given in [19], Appendix G, 
is best fit by r = 0.53. 

To fit the APM angular correlation function using a power spectrum for galaxies 
described by eq.(44) requires 0.15 ^ F < 0.3 [1] for rig = I and 0.2 ^ ris ^ 0.6 for F = 0.53 
[2, 74]. A recent estimate of F using power spectra from redshift surveys as well as from the 
APM data suggests F 0.25 fits best [75]. Fig.l. compares the COBE-normalized = 1, F = 
0.5 linear density power spectrum with an = 1, F = 0.25 and an rig = 0.6, F = 0.5 spectrum. 

To lower F into the 0.15 to 0.3 range one can: lower h; lower ilnr', or raise ilgr 
(= 1.69S1-^ with the canonical three massless neutrino species present). Low density CDM 
models in a spatially fiat universe (i.e. with A > 0) lower il„r to 1 — SIa- CDM models with 
decaying neutrinos raise Qgr [20, 76]: F 1.08Sl„rh(l -|- 0.96(myr(j/kev yr)^/^)"-'^/^ where 
is the neutrino mass and is its lifetime. Decaying neutrino models have the added feature 
of a bump in the power at subgalactic scales to ensure early galaxy formation, a consequence 
of the large effective il„r of the neutrinos before they decayed. 

Generally, more scales are needed to characterize the spectrum than just kneq, £-g-, the 
free-streaming scale for light neutrinos. In Hot/Cold hybrid models, there is a stable light neu- 
trino of mass rriu contributing a density Sl^ = 0.3(my/7.2 ev)(2h)~^, combining with the CDM 
and baryon densities to make a total il„r = 1. A F-shape is not a very accurate representation 
of the entire spectrum, dropping from about 0.5 for small A; to F ~ 0.22(51^/0.3)"^/^ over the 
band 0.04-2 h~"'^Mpc of relevance to Wgg(0) calculations [20, 77]. 

To lower rig, one can invoke one of the infiation models of § 2. utilizing a deceleration 
parameter q — («s -|- l)/2 or, for natural infiation, the curvature in InH away from the peak 

of the potential, ^ \nH/d(f)^ ^ (n,, - l)/2. 

Fig. 9. shows rip eff = dliiVp(k)/dliik for the models we have been discussing and 
compares it with the indications from observation, the 0.15 ^ F < 0.3 'box' extending from 
k~^ = 100h~^Mpc downward into the nonlinear region (indicated roughly by the intersection 
of the 7 = 1.8 line with 'box'.) Within the 'box', the preferred 0.25 value is also drawn. It can 
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be seen that the effective index function for the standard CDM model is well outside of the 
box. So also are the hot/cold hybrid models shown (with a 2.4 and 7.2 ev i>), although the 
Wgg{9) doesn't look too bad. The model with rig tilted to 0.5 falls within the 'box', as does 
a model with the lesser tilt, 0.7, but with Hq lowered as well, to 40. Of course, any rig = 1 
F-model with F ~ 0.25, such as one with a nonzero cosmological constant or a decaying v also 
falls within the box. The large circle shows the rough value for the slope indicated by the X- 
ray temperature distribution function of clusters and the smaller one, by the X-ray luminosity 
distribution function, although much can go wrong in the interpretation of the latter [78]. 

As we have seen, for inflation models we expect that there will be some tilt, so probably 
we can relax the F requirement. It is interesting to see what we would have to do to mimic a 
F = 0.2 shape with an initial spectral index variation if we were to assume a standard CDM 
model. The dotted ng(k) curve in Fig. 9. shows that radically-broken scale invariance is needed 
that changes the tilt from > 0.9 beyond ~ 200h-iMpc to ~ 0.5 at ~ 20h-iMpc. 



5.2. Relating the Cluster- ampHiude ag and the dmr Band-power 

Apart from the shape parameters for Vp(k), there is also an overall amplitude parame- 
ter, which we now take to be {C'f ^)^^^ = {Ci)^^^/(l + rig), eq.(31). Pre-COBE it was taken to 
be as, the rms (linear) mass density fluctuations on the scale of 8 h~"'^Mpc, which translates to 
the mass of a rich cluster, 1.2 x lO^^Q„r ((2h)~"'^ Mq: the number of rich clusters is extremely 
sensitive to the value of ag- Cluster X-ray data implies 0.6 ^ (Jg ^ 0.8 for CDM-like il„r = 1 
theories ([78], and references therein). We shall consider 0.7 as the target number and values 
below 0.5 as unacceptable. A higher value is better for SIa 7^ [1, 79]. 

The relation between ag and the dmr band-power is'* 

F-law: ^8 « 1-2 ,1 ^^-^'Imm 5 »n ° ""(2(r - 0.03)) e"^ «^(*-"') , (45) 
(1 + rtg)'-i^W ^ 

in \l/2 

hot/cold: ,78 « (Xn , (1 + 0-55(0^0.3)*/^)-* e-^-«^(i-"°) . (46) 

(1 -I- rtgyl^lQ--' 

This shows that for the observed band-power, the rig = 1 CDM model gives <Tg too 
high, 1.1, but this is a sensitive function of rig, dropping to 0.5 at rig = 0.82 with the standard 
gravity wave contribution or at rig = 0.7 if there is no gravity wave contribution. For the 
decaying neutrino model with rig = 1, to have <Tg > 0.5, we need F > 0.25 or m^Td < 7.5kevyr. 
For the hot/cold hybrid model, we need Sl^ 0.37 for <Tg = 0.7. See also ref.[80]. 

Another way to constrain the amplitude of the power spectrum is from the redshift 
of galaxy formation. We do not know this, but it cannot be too low or we would get too 
few z ~ 4 quasars and too little neutral gas compared with that inferred using the damped 
Lyman alpha systems seen in the spectra of quasars. This suggests 2 < (Tq^s ^ 5 or so, where 
(To 5 is the analogue of <Tg but at a Galactic mass scale rather than a cluster mass scale. This 
is the Achille's heel of hot/cold hybrid models [20, 80]. It also leads to serious constraints 
on rig for F = 0.5: in [2], we showed that a fairly conservative estimate of the redshift of 
galaxy formation was Zgj k. 1.3(To,5 — 1, and that (Tq^s ~ &.2(7ge~^^~"''^ for F = 0.5, leading 
to rig > 0.76 with gravity waves, rig > 0.63 without. With F < 0.5, the restrictions on the 
primordial spectral index from galaxy and cluster formation are even more severe. That is, 
only a little tilt is allowed. 

The free-streaming velocities of galaxies allow another estimate of the amplitude. For 
the F-laws, the large scale rms bulk flows over 40 and 60 h-*Mpc in km s"* are 

(7,(40 h-*Mpc)/,T8 = 317(2F)-°-^Xr '[1^0.572] e'"'('""') c./. 388 [1± 0.17] , (47) 

(7,(60 h-*Mpc)/,T8 = 252(2F)-°-^Xr '[1^0.572] e' ^('""') c./. 327 [1± 0.25] . (48) 



4 These fits were made in Ug with Y fixed at 0.5, and in Y with Ug fixed at 1, hence wih be rough when both 
vary signficantly from these standard values. 
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The comparisons are to the 'POTENT' estimates from the data [21]. The ±1(7 range is 
large because of the 'cosmic variance' expected in the measurement of a single bulk streaming 
velocity, but even so, the constraints in F — il„r — space are notable. For example, the 
F = 0.5 model needs rig > 0.89 with the typical gravity wave contribution and > 0.72 with 
none [2]; «,e,,in spite of the fact that tilted spectra have more large scale power, hence for 
a given <ts would have larger scale flows, the normalization to dmr implies just the opposite, 
that the index must be very nearly scale invariant. In [1], we used a different version of the 
velocity constraint with (formally) lower error bars: us ~ 1.0 ± 0.24 S1~^°'®®(2F)'^'^ for Ug = 1 
(derived by relating redshift survey results to streaming velocities [22], but with the simplifying 
assumption of a linear ampliflcation bias (Tg ^ for galaxies). This gives a (Ts incompatibity 
with the COBE-estimated value for = 1 A > models with F ~ 0.25. The peculiar 
velocity data relies on having spatially-independent and accurate distance indicators (e.^.jthe 
empirical Tully-Fisher relation between luminosity and rotation velocity in spiral galaxies). 
How seriously we take these constraints depends upon how reliable we think the indicators are 
- a subject of much debate. 

// the shape of the density power spectrum over the LSS hand is now considered to be 
known, then this restricts the range of inflation and dark matter models considerably, especially 
when combined with the COBE anisotropy level. Whether the solution is a simple variant 
on the CDM+mflation theme involving slight tilt, stable ev-mass neutrinos, decaying >kev- 
neutrinos, vacuum energy, low Hq, or some combination, is still open, but can be decided as 
the observations tighten, and, m particular, as the noise m the Ci figure subsides, revealing 
the details of the Doppler peaks, a very happy future for those of us who wish to peer into the 
mechanism by which structure was generated m the Universe. With the explosion of literature 
on such a wide-ranging topic, my referencing is far from complete - apologies. This work was 
supported by NSERC at Toronto and the Canadian Institute for Advanced Research. 
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Fig. 1. Cosmic waveband probes. The bands of cosmic fluctuation spectra probed by LSS 
observations are contrasted with the bands that current CMB experiments can probe. 
The (linear) density power spectrum for the standard = 1 CDM model, labelled 
r = 0.5, is contrasted with power spectra that fit the galaxy clustering data, one tilted 
(rij = 0.6, r = 0.5) and the other scale invariant with a modified shape parameter 
(rij = 1, r = 0.25). Biasing must raise the spectra up (uniformly?) to fit into the hatched 
Wgg range and nonlinearities must raise it at A; > 0.2hMpc~^ to (roughly) match the heavy 
solid (7 = 1.8) line. 

Fig. 2. (a) COBE-normalized temperature power spectra: for a standard = 1 CDM model 
with standard recombination, early reionization, a (dashed) tilted primordial spectrum 
with Hg = 0.95, with the gravity wave contribution shown, and a (dotted) high Hq model 
with A 7^ 0. Bandpowers with 10% error bars are shown for the tilted and untilted CDM 
models, (b) Filter functions for current experiments. 

Fig. 3. How a CDM model normalized to COBE varies with resolution. The contours begin 
at 109//A' in the half-degree smoothing cases, 54.5//A' in the no-smoothing case, 27.3//A' 
in the all-sky aitoff projection map. SR denotes standard recombination, NR denotes very 
early reionization, so there is no Doppler peak. The hills and valleys in the 5° SR map are 
naturally smooth: mapping them will give a direct probe of the physics of how the photon 
decoupling region at redshift ~ 1000 damped the primary signal. 

Fig. 4. Band-power estimates for the vintage Spring 1994 data, decribed in detail in the text. 
The lower panel is a closeup of the Doppler peak region. 

Fig. 5. (a) shows a 100° diameter map centred on the North Galactic Pole of the full 
dmr 53a-|-b data, while (b) shows it after the Weiner-filtering, assuming a n^T -|- 3 = 1 
Ci spectrum, (c) and (d) show the same for the 90a-|-b GHz data. The contours are 
±20n//A', n > 1 for (a,c), ±10n//A' for (b,d). Positive contours are heavier than negative 
ones. 

Fig. 6. (a,b,d) are the South Galactic Pole versions of Fig. 5. (c) shows the Weiner-filtered 
53a-|-b map with n^T -|- 3 = 2 is very similar in structure to (b) with n^T -|- 3 = 1. 

Fig. 7. Correlation functions for dmr (a-c) and firs (d-f) maps. (a,d) shows the correlation 
function for the data with k < (16)^ (closed triangles) and k > (16)^ (open circles and nearly 
zero, except within the beam). The x 's denote no cut, and for (a) the short horizontal lines 
denote k < (13)^. (b,e) are two simulated noiseless skies for n^T -|- 3 = 1 (open circles) 
and 2 (closed squares), with the same random number seeds. The k < (16)^ correlation is 
identical, while the k > (16)^ one has no correlation. (c,f) shows the correlation function 
for pure noise (x's: no cut; solid triangles: k < (16)^; open circles: k > (16)^). 

Fig. 8. Correlation functions for unfiltered maps (open circles, dotted error bars) are compared 
with those for the Weiner-filtered maps (open squares for n^T -|- 3 = 1) in (a)-(e). (a) and 
(c) also show riAT + 3 = 2 correlations (solid triangles). Maximum likelihood amplitudes for 
the band power were chosen in each case, (f) compares the rfmr Weiner-filtered correlations 
directly: 53 (open square), 90 (closed triangle), 31 (closed circle). The high noise maps are 
of course more filtered than the low. 
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Fig. 9. npeff(k) as a function of scale. The 0.15 
can get by adjusting cosmological parameters 
or by having 0.5 over the region of large 



^ r < 0.3 box is the target shape which one 
to ensure kjj]^^ ~ 15— 30 h~^Mpc for n,, = 1 
scale clustering observations for F = 0.5. 
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